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ABSTRACT 


HST/NICMOS images of the class I protostar TMR-1 (IRAS04361+2547) 
reveal a faint companion with 10.0" = 1400 AU projected separation. The 
central protostar is itself resolved as a close binary with 0.31" = 42 AU 
separation, surrounded by circumstellar reflection nebulosity. A long narrow 
hlament seems to connect the protobinary to the faint companion TMR-IC, 
suggesting a physical association. If the sources are physically related then we 
hypothesize that TMR-IC has been ejected by the protobinary. If TMR-IC has 
the same age and distance as the protobinary then current models indicate its 
flux is consistent with a young giant planet of several Jovian masses. 


Subject headings: binaries; general — circumstellar matter — infrared: stars 
planetary systems — stars: individual (TMR-1) — stars: formation 
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1. Introduction 


The past few years have seen the indirect detection by their gravitational effects of 
ronghly one dozen extrasolar Jupiter-mass planets around nearby stars ( [Marcy fc Butler 
|1998|) . Doppler surveys are primarily sensitive to giant planets within 3 AU of the central 
star, posing a challenge to theories which predict birthplaces in the 5-10 AU range. The 
conventional picture proposes a two-step formation process where a rocky planet core forms 
in a disk, followed by gas accretion over a period of 1 - 10 million years in a region outside 
5 AU radius ( |Lissauer 19951) . An alternative theory argues that giant planets form via 
gravitational instabilities in the disk on a time scale of thousands rather than millions of 
years ( Poss 199'^) . 

We have detected a low-luminosity object near the class I protostar TMR-1 whose 
flux is consistent with that of a giant protoplanet. If conhrmed, the protostar’s age of 
approximately 300,000 years places severe constraints on the time scale of giant planet 
formation. Furthermore the data show that the TMR-1 protostar is a binary system. Most 
stars are members of multiple star systems with separations ranging widely from sub-AU 
to thousands of AU; the overlap with circumstellar disk sizes has consequences for planet 
formation (Penest 1996 ; Pate fc Bonnell 1997|) . An important issue is whether binary stars 
provide a hospitable environment for the formation of substellar mass objects. 


2. TMR-1 protostars 

Previous observations establish TMR-1 as a typical class I protostar, similar in mass 


and luminosity to the Sun having ~ 0.5 Mq and 3.8 Lq, respectively, and which millimeter 
observations suggest has a low-mass disk (Pontemps et ah 1996|) . Observed NIR magnitudes 
are J = 16.1, H = 12.9, and K = 10.6 QTerebey et ah 19^ ; [Kenyon et ah 1993]) . Class 
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I protostars are surrounded by opaque envelopes of infalling gas and dust. Based on 
statistical and theoretical arguments typical ages of class I protostars are 100,000 to 300,000 
years (|Terebey et al. 1984| ; |Shu, Adams, fc Lizano 1987|) . NIR imaging and millimeter 
interferometry data of TMR-1 show a bipolar outflow which extends SE to NW (P?erebey et| 
|ah 19901 ; [Hogerheiide et al. 1998|) . Based on the arguments given in Chandler et al. (1996) 
TMR-1 is not viewed edge-on or pole-on, but at an intermediate 60°) inclination. 


3. HST NICMOS Images and Photometry 

The high spatial resolution (0.15" at 1.6 pm) HST/NICMOS images in Figure 1 resolve 
the TMR-1 protostar into two point sources which we call A and B. A is the northern 
component. At the Taurus cloud distance of 140 pc the 0.31" projected separation is 42 
AU, a fairly typical binary separation. The new data reveal TMR-1 to be a protobinary 
surrounded by gas and dust, viewed during the epoch of formation. 

Figure 1 displays extensive nebulosity, brightest near the protobinary. A long narrow 
filament extends in a gentle curve from near the protobinary to a third fainter point source, 
which we call C, located 10.0" southeast. The image provides strong visual evidence that 
object C appears associated with the protobinary by means of the filament. TMR-IC is 
detected at S/N = 50 in the F205W filter, as implied by the presence of the Airy diffraction 
ring. 

The image artefacts ([Casertano 1997|) are easily identihed in the original image 
orientation (-f-y axis at P.A. = 38° E of N). Artefacts arising from the bright protostars 
include: the ±45°telescope diffraction spikes; electronic ghost stars at ±128 pixels along 
the x,y axes; and two faint electronic ghost columns, one of which passes through the 
protostars, and another seen 128 pixels to the left. Finally, a residual coronographic spot 
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appears in the upper left hand quadrant. 


Extended nebulosity is common around protostars at near-infrared wavelengths. Stellar 
photons escape through the transparent polar regions created by the bipolar outflow; they 
delimit the r ~ 1 surface when they absorb or scatter in the dusty infall envelope and dusty 
circumstellar disk ([Whitney et al. 1997|) . The highly structured nebulosity around the 
TMR-1 protostars implies the density is inhomogeneous. The sharp contrast of the hlament 
above the background suggests the density is locally enhanced, while the illumination of the 
hlament at large (10") distances suggests a fairly clear line of sight back to the protostars. 


Table |I] gives positions and individual component huxes via PSF htting. Fluxes based 
on aperture photometry are given in Table ^ both to facilitate comparison with ground based 
measurements and to place the data on the STScI HST/NICMOS photometric system. 
Ground-based 2.2/im K-band IRTF images and also K' at Keck ( [Blake 1998|jHogerheijde et 
al. 1998[) conhrm the detection and approximate hux of object C. 


4. Local NIR Star Counts, Chance Background Object, and Extinction 

For the Taurus cloud K-band star counts give N{K' < K) = 0.041 x stars per 

square degree, which includes an extinction of Ak = 0.4 estimated from the same data 
( Beichman fc Jarrett 190^) . Assuming 18.5 for the K-band magnitude implies about one 
background star per NICMOS frame. 

A posteriori probability estimates are problematic. However we press on noting that 
the TMR-1 hlament is a unique structure in our ensemble of HST/NICMOS images. The 
chance that a random background star lies at the tip of the hlament is 2% if we assign a 
conservatively large 3" x 3" ehective search area. 


The scarcity of background stars is empirically conhrmed by the HST data, which show 
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fewer than expected background stars because of high extinction local to the protostars. 
Comparable S/N HST/NICMOS images for nine class I protostars in Taurus show one 
other secondary object (K = 18.7 mag), giving one or two possible background objects in 
nine helds. To match the large scale NIR star counts implies an average extinction over the 
20" NICMOS held of view of = 1 to 2 (Ay = 10 to 20) towards Class I protostars in 
Taurus. 


An alternate estimate for the extinction is set by values previously derived for the 
protostar, which range from 2.5 to 4 at K ([Terebey et ah 1990|; [Whitney et ah 1997|) . The 
extinction is likely smaller at 10" distance from the protostar, as is also suggested by the 
Table ^ hux ratios. Intrinsic NIR stellar colors are near zero because the spectral energy 
distribution of many stars is near the Rayleigh-Jeans limit. The observed highly reddened 
colors of protostars are therefore caused by extinction and scattering. The increasing hux 
ratio of object C to either protostar between 1.6 and 2.05 pm suggests less extinction 
toward object C than toward A and B. 


5. Luminosity and Temperature 


Models of giant planet and brown dwarfs imply they are hottest and brightest when 
young, as luminous as 0.01 Lq at one million years ( [Nelson, Rappaport, fc Joss 1993| ; 
Burrows et ah 19^) . The radii are near that of Jupiter’s, Rj = 7.1 x 10® cm, over a large 
mass and temperature range; young objects should be modestly (up to factor of three) 
larger. Models suggest ehective temperatures as great as 3000 K below one million years 
age. 


The object TMR-IC is clearly much fainter than the neighboring protostars; if located 


at the same distance as the Taurus cloud then the estimated bolometric luminosity is 
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approximately 10“^ to 10“^ Lq, within the giant planet to brown dwarf regime. To derive its 
luminosity the observed NIR fluxes were £t assuming for simplicity a black-body spectrum 
extincted by dust. Dust extinction parameters are from Draine (1998). Assumptions are 
140 pc distance, 1 Rj minimum radius for stellar and substellar objects, and Ay < 30 
{Kk < 3) extinction. The extinction cap is selected as the maximum compatible with NIR 
background source counts (Section^). 

General results from varying Tgj^, Ay, and radius are that the temperature is not well 
constrained as values ranging from 1200 K (Ay = 0) to 3000 K (Ay = 30) give acceptable 
hts. However the radius is reasonably constrained to be a few Rj at the assumed distance, 
and depending also on the maximum extinction at the high temperature end. Relaxing 
our assumptions, hotter background stars provide acceptable hts if higher extinctions are 
allowed. Foreground stars or low extinctions are ruled out; the approximate H - K color of 
1.5 is redder than the photospheres of known low luminosity stars. 


The broadband NICMOS filters give limited spectral information but allow us to 
exclude effective temperatures below approximately 1600 K since there is no evidence in 
TMR-IC for a strong methane dip near 1.8 microns (e.g. [AH^^rd fc Hauschildt 199^ ). 
Near-infrared spectra of cool objects (~ 2000 K) which show water at 1.9 pm are sufficiently 
featureless to be consistent with our photometry (e.g. Figure 7 of Pppenheimer et ah 1998|) . 
Better constraints on the extinction and effective temperature await low resolution spectra 
of TMR-IC. 


6. Mass 


Model evolutionary tracks for giant planets and brown dwarfs show the derived mass 
depends strongly on the age and luminosity (Figure 7 of Burrows et al. 19'T7 ). If TMR-IC 









has the same 300,000 year age assumed for the protostars then Ay is 8 - 20 and the mass 
is 2 - 5 Mj. If the age is ten million years, the same as older pre-main sequence stars in 
Taurus, the mass may be as high as 15 Mj. However, below one million years the models 
are sensitive to the initial conditions, as the thermal relaxation timescale is comparable to 
the planet’s age. More realistic models will depend on the planetary formation mechanism. 


7. Ejection Hypothesis 


If TMR-IC is a physical companion of the TMR-1 binary then models suggest it 
formed much closer to the protostars than its observed 1400 AU projected distance. We 
hypothesize TMR-IC has been ejected by the two protostars. Apart from some exceptions 
such as hierarchical systems, celestial dynamics hnds that 3-body stellar systems with 
comparable separations are unstable and tend to eject the lowest mass object (|Monaghan| 
1976 ). On dimensional grounds the characteristic velocity of ejection is (GM/R)°'^(1 -|- e). 


the velocity of periastron passage of the binary. Numerical studies show a large dispersion 
in ejection velocities (^tandish 1972|) . 


The separation of the protostars allows us to estimate a characteristic ejection velocity. 
The computation is only indicative given that the orbital parameters and inclination are 
poorly known. The observed projected separation of stars A and B is 42 AU; statistically 
binaries spend most time at the widest separations. For a typical binary eccentricity 
of e = 0.5 the separation varies by a factor of three. Including a modest deprojection 
correction, periastron passage may occur at 15 - 30 AU separation. The corresponding 
ejection velocity is 5 - 10 km s“^ for 1 Mq assumed total mass. The current distance of 10" 
then implies the time since ejection is about 1000 yr. 


Consider for the moment that the hlament marks the trail of object C. The hlament’s 
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shape is curved and appears consistent with the expected hyperbolic trajectory. However 
shear is likely important if the hlament lies within the differentially rotating infall envelope 
or disk. The assumption of Keplerian rotation is adequate to estimate the timescale 
( P?erebey et ah 1984|) . The period is 1000 yr at 100 AU radius, which implies signihcant 
wrapping can be expected on roughly 0.67" size scales. 


8. Filament 


Although the position angles of the hlament and outhow are similar, the hlament 
dihers from typical NIR outhow structures. Models of outhow cavities show conical shapes 
( [Whitney et ah 1997|) ; if the outhow cavity is limb-brightened it should have two symmetric 
horns with a sharp outer edge, whereas what is observed is one hlament whose sharp edge is 
on the wrong (southern) side given its curvature. NIR polarimetry data show the hlament 
is primarily scattered stellar light emanating from the protostars ([Lucas fc Roche 1997] ) 
which rules out an emission line jet. 


The hlament is projected against the outhow but the moderate source inclination 
implies the hlament could traverse either the outhow cavity or the dense infall envelope. 
One possible explanation is the hlament may be a material tail, such as for example a tidal 
tail formed by two colliding circumstellar disks ([Lin et al 1998[) . Ground-based data in 
support of a material tail show HCO’*' along the hlament, indicating the presence of dense 
gas ( [Hogerheijde et ah 1998[) . However HCO'*' can be ambiguous as a dense gas tracer since 
it often has enhanced abundance in molecular outhows. Alternatively the hlament may be 
an illumination channel, or light pipe, created when the protoplanet tunneled through the 
infall envelope. A drawback to the light-pipe explanation is that Bondi-Hoyle gravitational 
accretion implies a diameter which is too narrow to explain the observed hlament so some 
other mechanism must operate. 
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9. Isolated Planets 


We have proposed that TMR-IC is a substellar object which has been ejected by 
a binary protostar. There are two key experiments to test the idea that TMR-IC is an 
ejected protoplanet. Spectra will measure the extinction and effective temperature to better 
discriminate between stellar, brown dwarf, or planet origin. In several years proper motion 
measurements will detect TMR-lC’s motion on the sky. The predicted direction may be 
along, or in the case of a tidal tail, at an angle to the filament (|Lin et al 1998|). 


We outline one of the many possible mechanisms for planet ejection. Three-body 
numerical simulations suggest stable planetary orbits exist at radii approaching half the 
binary periastron separation ( [Benest 199(i|) . In other words there is a maximum stable 
radius for planet formation in a binary system. A substellar object that migrates or forms 
in the zone of marginal stability is subject to orbital resonance pumping. After repeated 
periastron passages the object gains sufficient energy to escape the system. This mechanism 
does not require a gaseous disk per se, and so may apply to pre-main sequence stars as well 
as protostars. 


The discovery of an ejected protoplanet is unexpected. However, given the prevalence 
of binary systems the process seems inevitable, and the question becomes how often. The 
idea that young planets should occasionally be ejected from their solar systems is rich in 
implications, both for our understanding of how planetary systems form, and in strategies 
for detecting isolated planets using current technology. 
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Fig. 1.— True-color NICMOS HST image of protostar TMR-1 with 3" = 420 AU scale. 
There is a long narrow filament extending from the protobinary to a faint point source 
which is the candidate protoplanet. The inset shows a magnihed view (0.3" scale) of the 
protobinary in the central region, complete with hrst Airy diffraction ring around each 
component. The held is suffused by bright nebulosity. The region is heavily reddened, bnt 
short wavelengths have been boosted to make the nebulosity appear white. Red (2.05 /im= 
F205W hlter); green (1.87 /im= F187W hlter); blue (1.60 /im= F160W hlter); North is up. 
Display uses log stretch. 



Table 1. TMR-1 HST NICMOS Positions® and Photometry^. 


Component 

Right Ascension 

Declination 

F160W 

F187W 

F205W 


(J2000) 

(J2000) 

1.60 ^m 

1.87 /cm 

2.05 /cm 

A 

04h39ml3.84s 

25°53'20.6" 

0 0+22 


250l?° 

B 

04h39ml3.83s 

25°53'20.4" 

q+17 

^-9 

23.5l^ 


C 

04h39ml4.14s 

25°53T1.8" 

1 1 7+0.06 
-0.04 


q qq+0.20 

O.OZ_o 20 


'^0.35" = Icr HST absolute position error. 

*’PSF amplitude fitting and Icr limits in Counts/s image units. Limits are 
instrumental for C, dominated by background structure for A and B. Absolute 
photometry errors of 3% should be added in quadrature. Applying no color correction, 
to convert units to /xJy multiply by 15.6, 34.4, and 12.4 for F160W, F187W, and 
F205W filters, respectively. 



Table 2. Flux of TMR-ICT 


R 

J 

F160W 

F187W 

F205W 

K 

L 

Units 



16. 

33. 

39. 



AiJy 

> 22.2'' 

> 21.2^= 

19.6 

18.5 

18.2 

17.9" 

> 13.2" 

mag 


^'NICMOS magnitudes are based on standard 0.5" apertures and the HST Vega 
system. Background structure dominates the 10% photometric uncertainty. 


5(7 limit from Jarrett et al. 1994 . 

April 1998 IRTF NSFCAM data. K band has 20% photometric uncertainty; 


3(7 limits elsewhere. 
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